Analysis of Scanning Tunneling Spectroscopy 
Experiments from First Principles: the Test 
Case of Ceo Adsorbed on Au(lll) ** 

Angel J. Perez- Jimenez,* Juan J. Palacios, Enrique Louis, 
Emilio SanFabian, Jose A. Verges 

February 6, 2008 



[* ] Dr. A. J. Perez- Jimenez, Prof. E. SanFabian 

Departamento de Quimica Ffsica, Universidad de Alicante 
San Vicente del Raspeig, Alicante E-03080 (Spain) 
Fax: (-1-34)965903537 
E-mail: AJ.Perez@ua.es 

Prof. J. J. Palacios, E. Louis 

Departamento de Ffsica Aplicada, Universidad de Alicante (Spain) 
Prof. J. A. Verges 

Instituto de Ciencia de Materiales de Madrid, CSIC (Spain) 

[** ] The authors wish to thank Drs. C. Rogero and J. I. Pascual for critically reading 
the manuscript and providing us with the experimental results shown in this paper. 
Financial support from Spanish DGICYT and the Generalitat Valenciana under 
Grants Nos. PB96-0085, GVOO-151-01 and GVOO-095-2 is gratefully acknowledged. 



1 



Many of the technological applications of fullerenes and their deriva- 
tives, either as lubricants, non-linear optical devices, superconductors or 
in molecular nanodevices depend on the use of thin films or the interac- 
tion with metal surfaces. This explains the tremendous effort that has been 
devoted by a number of researchers to understand and characterize Cgo ad- 
layers over a wide range of surfaces, usually by means of Scanning Tunneling 
Microscopy/Spectroscopy (STM/S) as well as with other spectroscopic tech- 
niques (See the work by Rogero et aM and references therein). Nevertheless 
STM/S data are sometimes difficult to interpret without aid from theoret- 
ical models. STM images and STS spectra are usually rationalized sim- 
plifying Bardeen's first-order formalismS by approximating the electronic 
structure of the tipil and that of the sample, as well as the tip-sample 
interaction.Bii On the other hand, single electron tunneling effects such 
as the Coulomb blockade and Coulomb staircase are accounted for by fits to 
simplified expressions derived from the golden-rule. Si Each of the afore- 
mentioned methods has its pros and cons, but none of them combines the 
following two highly desirable characteristics: First, to overcome Bardeen's 
approximation in favor of Landauer's formalism. This would allow a widen- 
ing of the analysis range beyond the tunneling regime down to tip-substrate 
contact. Second, to be fully ab initio, without resorting to ad hoc param- 
eters and approximations to describe the interaction between the adsorbed 
molecule, the substrate, and the tip. It is only in this way that a method can 
gain reliability and predictibility when reproducing STM/STS experiments. 



We propose here a method that brings these two major premises together 
and which has successfully been applied previously to the characterization 
of the conductance properties in molecular nanodevices of the type metallic 
electrode + molecule + metallic electrode. The procedure, which is 
termed Gaussian Embedded Cluster Method (GECM), computes Landauer's 
conductance formula through a divide-and-conquer scheme. This means that 
the electronic structure of the molecule and part of the electrodes is consid- 
ered in full depth from first principles, since it embraces the most important 
region, while the semi-infinite non-periodic structure of the metal is treated 
with an efficient recursive model. Herein we report its capabilities studying 
adsorbed Buckminster fuUerene on a gold surface, where recent and accurate 
STM/S experiments are available. A 

Nucleation and growth of Cgo monolayers and thin films on Au(lll) 
surfaces has been widely studied. 0El'0'll3Ei'01 Although binding between 
gold and Cqq is weaker than with other metals,000 it is far from neg- 
ligible; the adsorption energy is estimated around 40-60 kcal/mol. 

Bii In 

fact, adsorption of Ceo is able to lift the well-known 23 x -\/3 Au(lll) recon- 
struction, and photoemission studies of Ceo on polycrystalline Auil'ii and 
Au(lll) surfaces^ revealed energy shifts indicative of LUMO hybridization 
and charge transfer from Au to the adsorbed fuUerene molecule. At high cov- 
erages, closed-packed layers grow with the thermodynamically most stable 
adsorbate phase being a (2-\/3 x 2-\/3) R30" structure with a nearly perfect 
lattice matching in which all the molecules are in equivalent surface sites.© 
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Apart from this, another superstructure forms with crystallographic direc- 
tions matching those of the substrate, resulting in fullerene molecules sitting 
on different adsorption sites. 00 The proposed 11 x 11 Cgo coverage[El'0 for 
this superstructure has been recently discovered to be composed of a smaller 
2x2 grid by Rogero et alM In the following, we will focus on the character- 
istics of this superstructure as interpreted by Rogero et alM which serve as 
an excellent benchmark for our calculations to be compared with. 

The main features that characterize the adsorption of the above-mentioned 
Ceo adlayers can be drawn from Figure |I|, where we show the STS data 
obtained by Rogero et alM from Current Imaging Tunneling Spectroscopy 
(CITS) measurements. This technique measures the tunneling current at 
each point during a topographic scan for a range of bias voltages applied to 
the sample-tip system. The numerical differentiation of each curve provides 
the corresponding conductance profile. From these curves it is possible to 
construct an image of the conductance at a given bias for each point like 
the one appearing in the right panel of Figure ^ which was obtained at a 
bias of +O.6V.0 The curves plotted in the left panel of Figure |l| correspond 
to the conductance profiles measured at the points indicated with the circle 
and triangle. As deduced from Figure |l| the differences in brightness be- 
tween adsorbed Ceo molecules in the CITS image correspond to differences 
in conductance height of the peaks appearing in the conductance profile. 
On the other hand, the relatively weak C6o-Au(lll) binding reflects on the 
sharp form of the peaks, with the energetic position of the Ceo LUMO level 
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shifted towards the Fermi level by charge transfer from the substrate, with a 
HOMO-LUMO gap of 2.3 eV. The two curves appearing in the left panel of 
Figure |I| are interpretedS as corresponding to a 2x2 superstructure of the 
adsorbed layer that places the Cgo molecules on two different adsorption sites 
of the underlying substrate surface (see Figures |I|, |^ and p. According to this 
picture, Ceo molecules sit alternatively onto bridge and on-top sites, the dif- 
ferent interaction of the fuUerene with the adsorption site being responsible 
for the two type of spectra. Henceforth, when a Cgo molecule sits on top of 
a gold atom, where the interaction is supposed to be weaker, we end up with 
the bright peak, which is sharper than that corresponding to the molecule 
sitting on a bridge between two gold atoms, where the stronger interaction 
with the substrate would explain the increased width of the peak. The above 
site-specific properties are consistent with a larger displacement to the Fermi 
level of the peak corresponding to the bridge site with respect to that of the 
on-top one since the amount of charge transferred to the molecule is larger 
in the former than in the latter. 

The characteristics and interpretation outlined in the preceding para- 
graph have been confirmed by our calculations as can be seen after inspecting 
Figure |^ where we plot the conductance profile corresponding to the two ge- 
ometries that characterize the 2x2 Ceo superstructure. Since the interaction 
between the molecule and the surface is not very strong, the conductance 
peaks reveal the underlying positioning of the Ceo orbitals, from which we 
estimate the HOMO-LUMO gap to be 2.9 eV, in good accordance with the 
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experimental STS results. 10 The relative positioning and height of the con- 
ductance maxima between the on-top and bridge geometries is also repro- 
duced. The form of the peaks reflects, as suggested by Rogero et a/.,0 the 
different interaction of the molecules with the adsorption sites as deduced 
from the Potential Energy Scan (PES) shown in Figure |^. For each ad- 
sorption site we plot Density Functional Theory (DFT) calculations when a 
Ceo molecule approaches the Au(lll) surface towards the on-top and bridge 
sites, respectively, with either a six-member ring or a five-member ring fac- 
ing the surface. The same equilibrium distance (2.75A) is obtained for all 
the geometries which coincides with that found in Ceo-gold nanobridgesi^l 
and is consistent with Altman and Colton suggestion of no height differences 
between adsorbed Cgo molecules on Au(lll).tE10 As seen from Figure ^the 
adsorption energy is smaller when the molecule sits on top of a gold atom 
(16 kcal/mol) than when the fullerene binds to a bridge site (35 kcal/mol), 
which explains the narrower and higher form of the conductance maximum 
corresponding to the on-top geometry with respect to the bridge site. Once 
the aforementioned values are corrected by considering the lateral interaction 
of the Cgo monolayer £1, estimated from the Lennard- Jones potential to be 
of about 25 kcal/mol@, they lie within the experimental margins discussed 
above. The stronger Au-Cgo interaction when the molecule sits on the bridge 
site is consistent with a larger amount of charge transferred from the gold 
surface to the molecule: 0.8 electron, as compared to the value of 0.5 electron 
obtained for the on-top site. This gap accounts for the different alignment 
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of the Fermi level, which for the bridge geometry lies closer to the LUMO- 
derived orbitals, in complete agreement with Rogero et al. findings. The 
amount of transferred charge calculated by us is in accordance with recent 
results derived from photoemission spectraS, which give 0.8±0.2 electrons 
per fuUerene molecule adsorbed on Au(lll). Finally, we point out the fact 
that there exist minor differences in binding energy with respect to the sym- 
metry axis that points towards the surface. Actually, this coincides with the 
fact that no predominant molecular orientation had been found0 and is also 
indicative of the relative weak C6o-Au(lll) interaction, which allows a large 
degree of rotational freedom on the adsorbed fullerene. 

As mentioned before, our method is not only applicable to the tunneling 
regime but also at tip-sample contact distances. This is reflected in Figure ^ 
where we plot the conductance of the system for the two types of adsorption 
sites as the tip moves towards the adsorbed Ceo- The two maxima correspond 
to conductance channels coming from the first three Ceo LUMO-derived or- 
bitals, whose degeneracy has been partially removed due to the interaction 
with the gold surfaceBEl in two sets of two (broader peak) and one (sharper 
peak) resonances, respectively. These two peaks also appear at tunnel tip- 
surface distances in the experimental STS data of Rogero et al. (see Figure 

We can check that, as the tip comes in closer contact with the fullerene, 
their interaction alters the size, width and positioning of the peaks. The 
net result is an increase in the conductance to non- negligible values and a 
larger shift of the LUMO-derived Ceo orbitals towards the Fermi level due to 
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the larger amount of charge transferred from the tip atoms to the fullerene 
as they approach each other. This fact explains the slight closure of the 
HOMO-LUMO gap as the tip approaches the molecule, which is apparent 
after inspecting the insets of Figure |^. On the other hand, the relative shape 
of the peaks corresponding to the on-top and bridge geometries found at 
tunneling distances is maintained as the tip approaches the sample. This 
still reflects the differences in binding strength commented above. Figure ^ 
shows the value of the conductance maxima of the second peak vs. tip-surface 
distance. We can see the change of slope due to the different conductance 
regimes, with the exponential decay typical of tunneling appearing beyond 
13A, which agrees well with the experimental value reported by Joachim et 
alM of 13. 2A. The departure from a linear trend in the logarithmic repre- 
sentation of the conductance is related to the above-mentioned closure of the 
Ceo HOMO-LUMO gap. Structural deformation of the fullerene cage has 
not been considered in our calculations since we have focused on the change 
from tunneling to contact regimes, where the deformation of Ceo by the tip 
is negligible. 

In summary, we have shown that the so-called GECM method is able to 
accurately reproduce STS spectra. The example of Ceo adsorbed on Au(lll) 
represents just a good and difficult starting point, but the method looks 
promising as a valuable tool in the interpretation of STM and STS spectra 
up to the contact regime. 
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Methods 



Here we give a brief overview of our method; the interested reader is ref- 
erenced to our earher work. Our procedure aims at a first-principles 
implementation of Landauer's conductance formula: 

g = ^TT[TL&TRG'']. (1) 

In this equation, the operators Tl and Tn are built from the corresponding 
retarded, S^, and advanced S", self-energy operators of the left (L) and right 
(R) electrodes according to: 

^LiR) = ^{^L{R) - ^L(i?))> (2) 

while and G"- denote the retarded and advanced Green's functions opera- 
tors, respectively, of the whole system and where Tr represents the trace over 
the pertinent orbitals. Simple as it looks, Landauer's formalism poses a chal- 
lenging problem since the calculation of S involves the electronic structure 
of a semi-infinite non-periodic metallic electrode while G must deal with the 
intricacies of the metal-molecule interaction which, in turn, largely depends 
on the type and positions of the atoms involved. 

Our approach to the problem consists of a divide-and-conquer scheme. 
First, we perform a DFT calculation of the molecule including part of the 
leads with the desired geometry: see Figures |] and ^, where the two clusters 
representing the relevant part of the bridge and on-top geometries discussed 
above are depicted. We mention in passing that the DFT calculations were 
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performed by using the GAUSSIAN980 code with the B3LYP exchange- 
correlation functional© and the basis sets and Pseudopotentials of Chris- 
tiansen et a/.ii'il This is also the case for the DFT calculations that led 
to the PES of Figure ^ where clusters similar to those in Figures ^ and ^ 
were used, but with the five tip atoms removed. Additional clusters with 
the five-fold axis perpendicular to the surface were also used to derive the 
corresponding PES of Figure ^ but are not shown for convenience. 

Since the Hamiltonian of the cluster representing the nanobridge, H, is 
finite, its associated Green's functions are unsuitable for any current deter- 
mination. In order to transform this finite system into an effectively infinite 
one we must include the self-energy of the electrodes: 

[{e + i6)i-H-±'\eW{€) = i (3) 

where 

±^ = ±^^ + ±1. (4) 

The added self-energy is determined through a Bethe lattice tight-binding 
modelilii which is constrained to reproduce the electrode bulk density of 
states and to have the same Fermi energy as that of the system on which 
the DFT calculation was initially performed. The advantage of choosing a 
Bethe lattice resides in that the self-energies can be easily calculated through 
a well-known iteration procedure, il'il For each atom in the outer planes of 
the cluster we choose to add a branch of the Bethe lattice in the direction 
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of any missing bulk atom. Once the Green's functions have been calculated 
we proceed to self-consistency by recalculating the elements of the density 
matrix P^i according to: 



where ei;' is the Fermi level fixed by ensuring neutrality in the cluster, with S 
and G"" being the overlap and retarded Green's function matrices evaluated 
in a non-orthogonal atomic basis set. The new density matrix is used to 
recalculate the finite cluster Hamiltonian matrix, and the whole process is 
repeated until self-consistency is achieved. 

A final remark regarding the comparison between the experimental con- 
ductance profile and the one calculated by us must be made. Strictly speak- 
ing only the conductance at the Fermi level should be comparable with the 
experimental one, but the weak tip-Cgo interaction guarantees that the con- 
ductance profile will not vary strongly with the applied voltage within the 
bias range used in the STM experiments. This fact allows for a comparison 
between curves in the left panel of Figure |I] and those of Figure |^, although 
at contact tip-fuUerene distances the applied bias may produce non-trivial 
voltage drops. Extension of our method to cope with generic non-equilibrium 
situations is currently being implemented. 




mn 
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FIG. 1: Left panel: Normalized conductance derived from CITS measure- 
ments by Rogero et alM for the two types of Cgo adsorption sites on Au(lll). 
Right panel: CITS image by Rogero et alM at a bias of +0.6V, where the 
difference in brightness between on-top and bridge sites is clearly seen. The 
circle and triangle indicate the points where the curves shown in the left 
panel were obtained. The Fermi level is set to zero. The above pictures are 
reprints of Figs. 4d) and 4c) appearing in The Journal of Chemical Physics, 
Vol 116, No. 2, pp. 832-836. 

FIG. 2: Conductance spectra calculated with our method for the on-top 
and bridge adsorption sites in the region around the Fermi level (here set to 
zero). The inset shows the same data for the on-top geometry on a wider 
energy range, where the fuUerene gap can be clearly appreciated. 

FIG. 3: Cluster model used to represent the relevant part of the tip-Ceo- 
surface system. The fuUerene is oriented with an hexagon facing a bridge 
site of the Au(lll) surface. 

FIG. 4: Same as in Figure |^ for the on-top geometry. 

FIG. 5: DFT Potential Energy Scan of a Ceo molecule approaching the two 
types of Au(lll) surface sites (On-top and Bridge). Results with an hexagon 
and a pentagon of the fuUerene facing the surface are included for complete- 
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ness. The energy of the separated Ceo + surface is set to zero. 

FIG. 6: Top paneh Conductance profile from tunnehng to contact regime of 
a Cgo molecule adsorbed on top of a gold atom in a Au(lll) surface. Bot- 
tom panel: Same as above for the Cgo molecule adsorbed on a bridge site 
in a Au(lll) surface. Distance (in Angstrom) between gold surface and tip 
displayed in the legend. The Fermi level has been set to zero in both sets 
of curves. The inset in both panels plots the same data on a wider energy 
range at tip-surface distances of 13.5 and 14.5 A. 

FIG 7: Maximum conductance from the second peak of Figure |^ for both 
geometries. The tunneling regime can be distinguished by the exponential 
decay in the conductance beyond ISA. The inset shows the same data with 
a logarithmic scale for the conductance. 
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Figure 5: 
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Figure 6: 
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Figure 7: 




